Astronomy & Astrophysics manuscript no. sps paper' accepted 


© ESO 2009 


June 17, 2009 





SPS: A software simulator for the Herschel-SPIRE photometer 

B. Sibthorpe^ P. Chanial^, and M. J. Griffin^ 

' UK Astronomy Technology Centre, Royal Observatory Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ 
e-mail: sibOroe .ac.uk 

- Astrophysics Group, Blackett Laboratory, Imperial College, Prince Consort Road, London SW7 2AZ, UK 

e-mail: p.chanial@imperial.ac.uk 
' Cardiff School of Physics and Astronomy, Cardiff University, Queens Buildings, The Parade, Cardiff, CF24 3AA 

e-mail: matt . griff in®astro .cf.ac.uk 

Received XXXX XX, 2008; accepted XXXX XX, 2008 

ABSTRACT 

Aims. Instrument simulators are becoming ever more useful for planning and analysing large astronomy survey data. In this paper we 
present a simulator for the Herschel-SPIRE photometer. We describe the models it uses and the form of the input and output data. 
Methods. The SPIRE photometer simulator is a software package which uses theoretical models, along with flight model test data, to 
perform numerical simulations of the output time-lines from the instrument in operation on board the Herschel space observatory. 
Results. A description of the types of uses of the simulator are given, along with information on its past uses. These include example 
simulations performed in preparation for a high redshift galaxy survey, and a debris disc survey. These are presented as a demonstration 
of the sort of outputs the simulator is capable of producing. 
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1. Introduction 

Software simulators are now used in a wide range of astron- 
omy instrument projects. They have many uses, from guiding 
the instrument design, through observation planning, and finally 
helping to understand the astronomical data. This paper de- 
scribes the software package developed to simulate the Spectral 
and Photometric Imaging Receiver, SPIRE (Griffin et al.|2008bl l 
photometer instrument on board the European Space Agency's 
Herschel space observatory ( Pilbratt 2008 1, scheduled for launch 
in 2009. This software is designed to simulate the behavior of 
the instrument and spacecraft, and return realistic data in a for- 



mat compatible with the Herschel data pipeline ( |Giiffin et al 
|2008a| l. 

SPIRE is a dual instrument, comprising a three-band submil- 
limetre camera and an imaging Fourier transform spectrometer 
(FTS). In this paper we will only be considering the photome- 
ter; for information on the spectrometer simulator see 'Lindner 
|et al.| ( 2004| . The SPIRE photometer uses aiTays of hexagonally 
packed feedhorn-coupled bolometric detectors, and has a field 
of view (FoV) of 4 x 8 arcminutes. The FoV is observed si- 
multaneously in three bands centred approximately at 250, 350 
and 500 fim, giving diffraction limited full-width-half-maximum 
(FWHM) beam widths of approximately 18, 25 and 36 arcsec- 
onds respectively. There are 139, 88 and 43 detectors in the 
250 fim, 350 fim and 500 yum arrays respectively, along with 
two 'dark pixels' - bolometers positioned outside the instrument 
FoV - and two thermistors located on each array. 

Maximising the aperture efficiency of the feedhorns requires 
an aperture coiTesponding to an angle of 2A/D on the sky, 
where A is the wavelength and D is the telescope diameter 
Consequently, the detector beams have an angular separation of 
approximately twice the FWHM beam size on the sky. As a re- 
sult, specific observing patterns - either jiggling or scanning - 



must be employed to achieve Nyquist sampling of the sky bright- 
ness distribution (Griffin et al. 2002). Jiggling is achieved us- 
ing a dual axis internal beam steering mirror (BSM), while fully 
sampled scanned observations are obtained by scanning the tele- 
scope at the specific scanning angle of 42.4° with respect to short 
axis of the bolometer an^ay ( Sibthorpe et al. 2006 ; Waskett et al. 
|2007| |. 

The bolometer signal time-lines are filtered on-board by a 5- 
Hz low-pass filter. After multiplexing, an appropriate 4-bit offset 
is removed from each data time-line before the final gain stage, 
after which the data are digitized and telemetered to the ground. 
The offset removal allows the data to be sampled with 20 bit 
accuracy using a 16-bit ADC. The timelines are reconstructed 
to 20-bit accuracy on the ground using the telemtered offset and 
digitized data time-lines. The photometer on-board signal chain 
is described in detail by Griffin et al. (200 8a). 

In Section 2 of this paper we describe the simulation soft- 
ware itself, and its operation. Details of the format of the simula- 
tor input and output data are also given. Section 3 presents some 
examples of how the SPS has been used so far. This demon- 
strates the kind of outputs the SPS can produce, and the types 
of investigations that it can be used to perform. The usage per- 
missions and availability of the SPS are explained in Section 3, 
and Section 4 describes the future development plans for the SPS 
post-launch. Finally, Section 6 reviews the main conclusions of 
this paper. 



2. The SPIRE photometer simulator 

The SPIRE photometer simulator (SPS) is a software package 
which simulates the operation of the SPIRE photometer instru- 
ment and associated spacecraft functions. The data are output as 
time-lines in engineering units, which corresponds to the input 
data level (Level-0) for the Herschel data pipeline. 
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The SPS allows users to input their planned observations, 
along with representative maps of the sky (either synthetic or 
derived from existing data), and be returned data which can be 
reduced and analysed using the official Herschel data pipeline, 
or any other suitable data reduction software the user wishes to 
use. 

It is important to note that the SPS should not be used as 
an observation planning tool as it does not take into account 
instrument and spacecraft overheads. All observation planning 
should be carried out using the Herschel observation planning 
tool HSpot ( Herschel-Spot Users' Guide 2007). The simulated 
data are intended to be a faithful representation of what in-flight 
data will be like, but while every effort has been made to make 
the data realistic, it is inevitable that as yet unknown instrument 
systematics will be present in the post-launch data which are not 
currently present in the simulated data. 

The SPS is controlled via a graphical user interface (GUI). 
Observations are set up using a window which contains the same 
parameters as those in HSpot. A simulated observation can then 
be performed using any of the SPIRE astronomical observing 
modes, including parallel mode (only SPIRE output is returned 
by a parallel mode simulation). 

2.1. Simulator architecture and operation 

The SPS is a modular software package written in the Interactive 
Data Language (IDL), with each module being controlled by a 
core program. This design was chosen to allow various modules 
to be developed independently based on a defined set of module 
input and output parameters. With the exception of the core pro- 
gram, each module contains four different phases of operation: 
an input, initialisation, temporal, and afterrun phase. The core 
module contains only the input phase operation. 

The input phase loads all of the user defined inputs and stores 
them in the appropriate data structures in preparation for the 
simulation run. The initialisation phase uses these input data 
to compute a variety of standard parameters, including optical 
transmissions, detector array layout etc. These are parameters 
which will not change throughout the simulations, and thus can 
be computed at the start of the simulation. Parameters whose 
values might change as a result of feedback from other system 
components are computed in the step phase. Most commonly 
these will be sets of parameters which are mutually dependent. 
In most cases the majority of a module's operation can be per- 
formed in the initalisation phase, with the step phase allowing 
for flexibility in the system for future development. Finally, the 
afterrun phase is used to compute any parameters which require 
entire data time-lines to be complete before their operation can 
be performed. For example, the 1// noise in the SPS is com- 
puted using Fourier methods, and cannot be computed in a step 
by step method. Likewise the simulated on-board filtering of the 
data must be done on entire data time-lines. 

The processing sequence of the SPS is shown in Figure[T]as a 
flow diagram, with each phase of operation being displayed. This 
diagram also shows how multiple simulation runs are performed. 

2.2. Modules 

The SPS contains 1 1 modules, with each module representing 
the operation of a specific system component. For example, the 
detectors module contains all code associated with simulating 
the detectors, including the generation of noise contributions 
arising from the detector itself, and the associated load resistor 
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Fig. 1. Block diagram showing the operational architecture im- 
plemented in the SPS. 



The modules are named according to the subsystem or function 
that they simulate and are described in detail below. 

2.2.1. Core 

This module does not represent any physical part of the 
Herschel-SPIRE system. Instead it manages and runs all oper- 
ations specific to operating the simulation software, such as call- 
ing modules in the correct order, performing the required num- 
ber of simulations, and setting the simulation output file names. 
It also controls the time step resolution used within the simula- 
tor, which directly relates to the simulation run time. 

The nominal time step used in the SPS is 21 ms, which 
equates to 48 Hz, approximately 3 times the nominal instrument 
detector sampling rate. The internal SPS clock must run faster 
than the instrument sampling rate to account for short time-scale 
transient systematics, such as high frequency noise fluctuations. 
The 3:1 ratio of instrument to internal clock sampling rate was 
found to give provide an accurate simulated output, with no im- 
provement found when operating at higher ratios. 

2.2.2. Sky 

Like the core, this module does not correspond to a particular 
instrument system. Instead it represents the astronomical sky 
viewable by the telescope. The module reads in three user de- 
fined input sky files, one for each waveband in a Flexible Image 
Transport System (FITS) file format. The units of the input maps 
are Jy/pixel, and the pixels must have the same integer arcsec- 
ond size for all three bands. It is recommended that a pixel size 
of 2 arcseconds is used (~1/10 of the 250 fim beam). Pixel sizes 
larger than this can result in digitisation artefacts from the in- 
put map being present in the output data. The source spectrum 
in each pixel is treated uniform across the SPIRE band, there- 
fore any spectral slope information must be included in the input 
model. This input will, however, be multiplied with the SPIRE 
filter response function. 
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The astrometric data associated with the input images are 
also read and stored for use in the observatory function module. 
The astrometic information contained within the FITS header 
must therefore be correct and contain the CDELT, CRPIX, 
CTYPE, CRVAL, CROT and EQUINOX keywords. 

2.2.3. Optics 

This module represents the main optical properties and param- 
eters of the telescope and the SPIRE photometer (including the 
filters), and the positional mapping of the detectors on the sky. 
These parameters are used to derive the optical transmission and 
emissivity profiles for each optical element, as well as the final 
instrument spectral response function in each band. This spec- 
tral response function is then applied to the input sky model as- 
suming a flat source spectrum across the SPIRE bands. The true 
source spectrum must be taken in to account in the generation of 



the input sky model (see Section 2.2.2 1 



The telescope beam is also taken into account in this module, 
and is used to produced a beam convolved version of the user's 
input map. A single indentical beam model is used for each de- 
tector within a single array, and three beam models are used in 
total, one for each band. The beams used are based on optical 
models and include the telescope side-lobes. The beam model is 
contained within a library file and loaded from memory when re- 
quired. Each beam is symmetric and is contained within a square 
grid of size -4x4 FWHM. 

2.2.4. Observatory Function 

The observatory function module specifies the photometer ob- 
serving mode to be simulated in terms of the appropriate ob- 
serving function and its associated parameters. Observing pa- 
rameters are supplied via either a simple HSpot style interface, 
or an advanced interface option. The HSpot option is the nom- 
inal input method, and provides input parameters akin to those 
found in the HSpot software. This does not mean however that 
the software can read in HSpot output files, rather the input pa- 
rameters are defined in the same way as those in HSpot. The ad- 
vanced screen allows the user to select different options in order 
to investigate the effect of different observing strategies which 
are not normally permitted by Herschel. 

The user input information supplied is then used to gener- 
ate both the telescope pointing information, including pointing 
errors, as well as the input parameters for the BSM module. 



a user specified noise spectral density, knee frequency, and spec- 
tral index (a). 

The cooler fluctuation time-line is low-pass filtered assum- 
ing an RC filter profile. This represents the suppression of high 
frequency thermal variations by the various thermal linkages be- 
tween the cooler tip and the array, and the thermal mass of the 
detector arrays themselves. A different user defined time con- 
stant is used for each array, representing the different distances 
and linkages between the cooler tip and each specific array. This 
results in the three arrays having a similar, although not identical 
thermal fluctuation time-line. A demonstration of the impact of 
thermal drifts on the output timeline can be seen in Figure [2]; a), 
along with the reference astronomical power timeline. 

At present there are no thermal gradients within a bolometer 
array, meaning that all bolometers vary in temperature simulta- 
neously. These thermal variations represent the dominant source 
of common mode 1/f noise in SPIRE, and as with all systemat- 
ics simulated, can be turned on and off depending on the users 
needs. 



2.2.7. Astronomical Power 

This module uses both the telescope and BSM pointing infor- 
mation to generate absorbed power time-lines for each detector. 
The observing strategy is projected onto the beam convolved in- 
put sky, and the flux from each pointing is registered. Additional 
telescope efficiency factors and filter transmission profiles are 
then used to compute the astronomical power falling on a given 
bolometer during each simulation time step. A sample output 
timeline is given by the dashed line in Figure |2] 



2.2.8. Background Power 

The background power falUng on each bolometer throughout 
a simulation is generated in this module. Each element in the 
optical chain is assumed to emit as a blackbody, with a wave- 
length dependent emissivity of (1 - t), where t is the optical 
transmission. The primary miiTor is modelled with a user de- 
fined thermal gradient (with temperature decreasing with dis- 
tance from the spacecraft's sun shield). The background power 
of the Herschel/SPIRE system is expected to be stable on the 
period of a single observation, therefore no background power 
variations are implemented. This will be reassessed once post- 
launch data are available. 



2.2.5. BSM 



2.2.9. Detectors 



This module simulates the movement of the BSM. It generates 
pointing time-lines in terms of the fixed array spacecraft axes. 
These time-lines are generated originally in arcseconds, and then 
converted to analog-to-digital units (ADU) via a look-up table. 
Both outputs are provided to the user. 

2.2.6. Tliermal 

This module simulates all information pertaining to the tempera- 
tures of the instrument and the telescope, and their temporal fluc- 
tuations. However, since the instrument system, and telescope 
can be considered to be stable on time scales of a single obser- 
vation, the model used in the SPS assumes that the only source 
of thermal fluctuations is the ^He cooler. The cooler fluctuations 
are simulated using a simple l/f" noise time-line generator, with 



This module produces an output voltage time-line for each 
detector channel based on inputs from the Astronomical and 
Background power modules. The bolometer thermal model pre- 
sented by |Sudiwala et al.| ( |T992| and |Woodcraft et al.| ( [2008| l is 
used to calculate the correct voltage output based on the power 
absorbed by the detector. The thermal model provides an accu- 
rate characterisation of non-linear bolometer response to high 
levels of incident power, and detector output fluctuations due to 
variations in the bath temperature. Each detector is modelled in- 
dependently using its own set of bolometer parameters measured 
during instrument testing. 

In addition to the detector voltage time-lines, the theoretical 
photon, Johnson, phonon and amplifier noise contributions are 
all calculated here. These data are then used to generate unique 
white and/or l/f noise time-lines, which are in turn summed 
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Fig. 2. Sample data timeline for the central bolometer in the 
250 fim array under various noise conditions: (a) contains only 
thermal drifts (correlated 1/f noise) with no other noise effects; 
(b) shows the timeline with uncorrected l/f and white noise ap- 
pUed; (c) contains both of the previous two outputs combined, 
and (d) shows the timeline from (c) once it has passed through 
the read-out electronics, i.e. low-pass filtered, sampled, and digi- 
tised. All signals have been calibrated to Jy/beam from their na- 
tive unit to allow direct comparison for direct comparison of the 
different noise effect. The original noiseless astronomical power 
timeline is also given in each panel by the dashed line. 



with the detector output to produce the final detector voltage 
time-line (see Figure|2|b)). 

2.2.10. Sampling 

The sampling module simulates the processing of the output 
detector time-lines by the on-board read-out electronics. Each 
detector time-line is filtered by the on-board 5-Hz low-pass fil- 
ter, digitized, and output at the requested sampling frequency. A 
voltage offset is also removed from the time-lines, and output as 
a separate single value for each detector. 

Sample input and output time-lines to this module are given 
in Figures |2|c) and (d) respectively. The decreased noise level 
resulting from the low-pass filtering, along with the digitisation 
of the signal are evident in these panels. 



2.2.11. Housekeeping 

The housekeeping module outputs various instrument monitor- 
ing information. These data include flags which identify the ob- 
serving mode being used, and the building blocks used to make 
up that mode. 



2.3. Output Data 

The SPS output is a standard FITS file with 11 extensions. 
Each extension contains the data structure from one of the SPS 
modules. All input parameters are contained within the output 
file, however for file size reasons not all derived time-lines are 
stored as standard. Nominally, only the time-lines returned by 
the physical system are output, i.e. telescope pointing and digi- 
tised bolometer outputs. It is possible, however, to choose to out- 
put the astronomical power, background power, noiseless volt- 
age, and noisy voltage time-lines if requested. These time-lines 
are sampled at a higher frequency, and hence can significantly 
increase the output file size. 

The standard FITS output can be loaded into the official 
Herschel pipeline processing software via a separate piece of 
converter software, whereupon it can then be reduced in the 
same way as real SPIRE data. Alternatively the standard output 
can be converted from digitized units (level-0) to calibrated data 
(level- 1) time-lines using a conversion routine supplied with the 
SPS. This new file can then either be loaded into the Herschel 
pipeline, or read in using another piece of custom software. 



2.4. System requirements, performance and limitations 

The SPS will run on any modern computer (CPU ~2 GHz, RAM 
~2 GB) and requires IDL V6.2 or higher Using a typical com- 
puter the SPS will run approximately 4 times faster than real 
time, meaning it would take 1 hour to simulate a 4 hour obser- 
vation. Typically a simulated one square degree scan map obser- 
vation will require ~635 MB of memory during the simulation, 
and generate a ~44 MB output file. The memory required during 
a simulation, and the size of the output file are both dependent on 
the size of the input file, and the selection of output parameters. 
In an optimal configuration for specific investigations the sys- 
tem has been shown to run up to 10 times faster than real time. 
Simulated observations can be performed for the full 18 hours 
available in a single Herschel astronomical observing request. 

A lack of detail in the spacecraft pointing model means that 
SPS cannot be used to calculate observing times. Many of the 
required overheads such as slewing and pointing calibrations are 
not included in the SPS. The pointing model used has been de- 
signed to provide realistic data for on-source regions only. 

The SPS contains all of the major known instrument sys- 
tematics, and represents the best current estimate of the SPIRE 
photometer's in-flight performance. 



3. Uses of the SPS 

The SPS has been used for a variety of purposes throughout its 
development. It was used to optimise and characterise the SPIRE 
photometer astronomical observing modes, and identify the best 
observing strategy possible with SPIRE ( Sibthorpe et al.|[2006 



IWaskett et al":][2007| [Sibthorpe et "aD '2008a'). It has also pro 
vided input data for a variety of map-making algorithms during 
a map-making selection exercise. The different algorithms were 
compared, and the best chosen as the standard map-making al- 
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gorithm for use in the Herschel pipeline (Chanial, submitted to 
PASP) 

In addition to mission planning, the SPS is currently be- 
ing used in the preparation and planning of various Herschel 
key programmes, both in guaranteed and open time. These pro- 
grammes span a very wide range of targets, from high redshift 
galaxy surveys through to observations of diffuse galactic struc- 
ture and star forming regions. Investigation of saturation in sur- 
veys containing high dynamic range, uniformity of noise in out- 
put maps, and development and testing of specific algorithms for 
use with survey data are all examples of recent uses of the SPS. 

An example output produced during planning for one of the 
high redshift galaxy surveys is presented in Figure [3] The input 
to this simulation is a map derived from galaxy catalogues which 
were scaled to the appropriate flux densities. Three input images 
were created, one for each band, but only the 250 fim band im- 
age is shown here. Two orthogonal simulated 'large map' scans 
were performed with 1// noise included. The time-line data from 
these two observations were then made into a naive map and co- 
added to produce the final image seen in Figure |3^. This type 
of image can be used to investigate the kind of problems this 
type of survey might have with source extraction, P(D) analysis, 
and various other data analysis techniques. A coverage map, also 
output by the map-maker (Figure [S]?), provides additional infor- 
mation on the uniformity of integration time, and hence depth 
across the map. With these two maps the noise statistics of the 
data can be investigated, and any biases characterised. With a 
known input sky model, methods which attempt to remove the 
various sources of noise can also be tested and compared. 

An example of a 'small map' observation is given in 
Figure|4] Here the observation of a debris disc, whose brightness 
distribution is based on that measured for the Epsilon Eridani 
disc ( [Greaves et al.|[T998] l, has been simulated without instru- 
mental noise. The input sky model is shown in Figure]?^, and 
consists of a simple torus shaped object placed on top of a real- 
istic extra-galactic background. The goal of this simulation was 
to provide a way to study the influence of extra-galactic confu- 
sion on the derived properties and structure of the debris disc 
when performing a chopped SPIRE observation. This investi- 
gation found that a 'small map' observation results in a ~20% 
higher extra-galactic confusion noise level than the same obser- 
vation performed in 'large map' mode ( Sibthorpe et al.|2008b i. 
In addition, it was found that the structure of the background 
galaxies in the 'small map' simulation was far more uncertain as 
a result of their emission being mixed in the output with that of 
sources in the off-source field. Consequently, methods to iden- 
tify and remove the contaminating sources are significantly more 
difficult to employ. 

The information from these types of simulations allows for 
a better informed choice of observing mode. A mode can be se- 
lected which fits the type of source, and confusion noise environ- 
ment. Simulations can also be used to develop, test, and charac- 
terise a series of analysis routines so that they are ready once the 
data are taken. 
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Fig. 3. Simulated 'large map' observation of a high redshift 
extra-galactic field, (a) This image is from the 250 /urn band, and 
shows an image of size 45 arcmin x 45 arcmin. A naive map- 
maker was used to construct the image from the output time- 
line data, (b) coverage map showing uniformity of the observing 
time. 



the models used within the SPS. This site will also contain in- 
formation on updates to the SPS as they are released'. 



4. Usage and availability 

The SPS software is publicly available from the SPS website 
(www.roe.ac.uk/~sib/sps), along with the required conver- 
sion software which enables the SPS output to be used with 
HIPE. All relevant information on how to obtain and use the SPS 
can be found at on the SPS website, along with further details on 



' The SPS may be used free of charge only for non-commercial re- 
search purposes. If you make modifications to this software, you must 
clearly mark the software as having been changed and you must also re- 
tain the copyright and disclaimer In the interests of avoiding the prop- 
agation of modified versions of the SPS, we would request that you do 
not edit or redistribute this software without first contacting the author, 
or a member of the SPS team. All uses of this software must acknowl- 
edge it's use, and reference this paper. 
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Fig. 4. Noiseless simulated 'small map' observation of a debris 
disc - (a) input sky model, (b) output map. 
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5. Future development 

Following the launch of Herschel a new version of the SPS will 
be released. This version will be updated to include in-flight per- 
formance data, and will aim to replicate the true characteristics 
of the operating system. Feedback from users is also expected to 
guide the development of additional features in the SPS. 



6. Conclusions 

In this paper we have described the SPIRE photometer simulator 
software and given examples of its capabilities and outputs. The 
SPS mimics the output of the SPIRE photometer such that it can 
be loaded into the Herschel pipeline environment, and reduced 
as if it were real data. The software represents the primary in- 
strumental characteristics, and allows the user to investigate how 
these characteristics might influence the quality of astronomical 
data for individual science cases. It also provides an opportunity 
to prepare, and characterise, science case specific data analysis 
routines. 

We have presented several examples of preparation work car- 
ried out with the SPS, and provided information on how new 
users can access and use the SPS in preparation for their own 
Herschel observations. The SPIRE photometer simulator is a 
versatile tool for the investigation of a wide range of instru- 
mental influences on astronomical data obtained with Herschel- 
SPIRE. 
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